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Precision measurements, now and at a future linear electron-positron collider (ILC), can 
provide indirect information about the possible scale of supersymmetry. We illustrate the 
present-day and possible future ILC sensitivities within the constrained minimal supersym- 
O \ metric extension of the Standard Model (CMSSM), in which there are three independent soft 



supersymmetry-breaking parameters mi/2, rn^ and Aq. We analyze the present and future 
p • sensitivities separately for Mw, sin^6'efr, {g — 2)^, BR(6 — > 37), BR(i?5 — > yU^/i"), Mh and 
^ ' Higgs branching ratios. We display the observables as functions of mi/2, fixing mo so as to 
obtain the cold dark matter density allowed by WMAP and other cosmological data for spe- 
^ . cific values of Aq, tan/5 and /i > 0. In a second step, we investigate the combined sensitivity 
■ of the currently available precision observables. My/, sin^6'efr, {g — 2)^ and BR(6 S7), by 
performing a analysis. The current data are in very good agreement with the CMSSM pre- 
diction for tan/? = 10, with a clear preference for relatively small values of mi/2 ~ 300 GeV. 
In this case, there would be good prospects for observing supersymmetry directly at both 
the LHC and the ILC, and some chance already at the Tevatron collider. For tan (3 = 50, 
the quality of the fit is worse, and somewhat larger mi/2 values are favoured. With the 
prospective ILC accuracies the sensitivity to indirect effects of supersymmetry greatly im- 
proves. This may provide indirect access to supersymmetry even at scales beyond the direct 
reach of the LHC or the ILC. 
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1 Introduction 



Measurements at low energies may provide interesting indirect information about the masses 
of particles that are too heavy to be produced directly. A prime example is the use of precision 
electroweak data from LEP, the SLC, the Tevatron and elsewhere to predict (successfully) the 
mass of the top quark and to provide an indication of the possible mass of the hypothetical 
Higgs boson [1]. Predicting the masses of supersymmetric particles is much more difficult 
than for the top quark or even the Higgs boson, because the renormalizability of the Standard 
Model and the decoupling theorem imply that many low-energy observables are insensitive to 
heavy sparticles. Nevertheless, present data on observables such as Myy, sin^ 6'eff, — and 
BR(6 — s> s'-f) already provide interesting information on the scale of supersymmetry (SUSY), 
as we discuss in this paper, and have a great potential in view of prospective improvements 
of experimental and theoretical accuracies. 

In the future, a linear e^e~ collider (ILC) will be the best available tool for making 
many precision measurements [2]. It is important to understand what information ILC 
measurements may provide about supersymmetry, both for the part of the spectrum directly 
accessible at the LHC or the ILC and for sparticles that would be too heavy to be produced 
directly. Comparing the indirect indications with the direct measurements would be an 
important consistency check on the theoretical framework of supersymmetry. 

Improved and more complete calculations of the supersymmetric contributions to a num- 
ber of low-energy observables such as M]y and sin^ ^efr have recently become available (see 
the discussion in Sect. El below). These, combined with estimates of the experimental accu- 
racies attainable at the ILC and future theoretical uncertainties from unknown higher-order 
corrections, make now an opportune moment to assess the likely sensitivities of ILC mea- 
surements. 

There have been many previous studies of the sensitivity of low-energy observables to the 
scale of supersymmetry, including, for example, the precision electroweak observables [3-9]. 
Such analyses are bedevilled by the large dimensionality of even the minimal supersymmet- 
ric extension of the Standard Model (MSSM), once supersymmetry-breaking parameters are 
taken into account. For this reason, simplifying assumptions that may be more or less well 
motivated are often made, so as to reduce the parameter space to a manageable dimension- 
ality. Following many previous studies, we work here in the framework of the constrained 
MSSM (CMSSM), in which the soft supersymmetry-breaking scalar and gaugino masses are 
each assumed to be equal at some GUT input scale. In this case, the new independent MSSM 
parameters are just four in number: the universal gaugino mass mi/2, the scalar mass mo. 
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the trilinear soft supersymmetry-breaking parameter Aq, and the ratio tan/3 of Higgs vac- 
uum expectation values. The pseudoscalar Higgs mass Ma and the magnitude of the Higgs 
mixing parameter /i can be determined by using the electroweak vacuum conditions, leaving 
the sign of // as a residual ambiguity. 

The non-discoveries of supersymmetric particles and the Higgs boson at LEP and other 
present-day colliders impose significant lower bounds on m.1/2 and mo- An important fur- 
ther constraint is provided by the density of dark matter in the Universe, which is tightly 
constrained by WMAP and other astrophysical and cosmological data [10]. These have the 
effect within the CMSSM, assuming that the dark matter consists largely of neutrahnos [11], 
of restricting m.Q to very narrow allowed strips for any specific choice of Aq, tan/3 and the 
sign of II [12,13]. Thus, the dimensionality of the supersymmetric parameter space is further 
reduced, and one may explore supersymmetric phenomenology along these 'WMAP strips', 
as has already been done for the direct detection of supersymmetric particles at the LHC and 
linear coUiders of varying energies [14-19]. A full likehhood analysis of the CMSSM planes 
incorporating uncertainties in the cosmological relic density was performed in Ref . [20] . The 
principal aim of this paper is to extend this analysis to indirect effects of supersymmetry. 

We consider the following observables: the W boson mass, Mw, the effective weak mix- 
ing angle at the Z boson resonance, sin^ ^eflf, the anomalous magnetic moment of the muon, 
{g — 2)^ and the rare b decays BR(6 — > 57) and BR(Ss — > ji^ji"), as well as the mass of the 
lightest CV-even Higgs boson, M^, and the Higgs branching ratios BR(/i — > 66)/BR(/i — > 
1^1^*). We first analyze the sensitivity of each observable to indirect effects of supersym- 
metry, taking into account the present and prospective future experimental and theoretical 
uncertainties. We then investigate the combined sensitivity of those observables for which 
experimental determinations exist at present, i.e., M^, sin^ ^eff, {g — and BR(6 57). 
We perform analyses both for fixed values of Aq and for scans in the (mi/2, ^0) plane for 
tan /3 = 10 and 50 with // > 0. We find a remarkably high sensitivity of the current data for 
the electroweak precision observables to the scale of supersymmetry. In the case tan /3 = 10, 
we find a preference for moderate values of mi/2 ~ 300 GeV, in which case sparticles should 
be observable at both the LHC and the ILC. In the case tan /3 = 50, the global fit is not 
so good, and low values of mi/2 are not so strongly preferred. In order to investigate the 
possible future sensitivities we study the combined effect of all the above observables (except 
BR(5s — > jjL^jjr), which is discussed separately). For this purpose we choose certain values 
of (mi/2, Aq) as assumed future 'best-fit' values (corresponding to the central values of the 
observables) and investigate the indirect constraints arising from the precision observables 
for prospective experimental and theoretical uncertainties. 



2 



In Section 2 of the paper we specify the WMAP strips and discuss their dependences 
on ^0 and the top-quark mass. We discuss in Section 3 the present and future sensitivities 
of the different precision observables to the scale of supersymmetry, represented by mi/2 as 
one moves along different WMAP strips. In Section 4 we analyze the combined sensitivity 
of the precision observables for the present situation, and Section 5 presents the prospec- 
tive combined sensitivity assuming the accuracies expected to become available at the ILC 
with its GigaZ option. Finally, Section 6 gives our conclusions. In most of the scenarios 
studied, even if it does not produce sparticles directly, the ILC will check the consistency of 
the CMSSM at the loop level and thereby provide valuable extra information beyond that 
obtainable with the LHC. 

2 Supersymmetric dark matter and WMAP strips 

It is well known that the lightest supersymmetric particle (LSP) is an excellent candidate 
for cold dark matter (CDM) [11], with a density that falls naturally within the range 
0.094 < flcBuh^ < 0.129 favoured by a joint analysis of WMAP and other astrophysical 
and cosmological data [10]. Assuming that the cold dark matter is composed predominantly 
of LSPs, the uncertainty in the determination of JIcdm^^ effectively reduces by one the 
dimensionahty of the MSSM parameter space. Specifically, if one assumes that the soft 
supersymmetry-breaking gaugino masses mi/2 and scalar masses mo are universal at some 
GUT input scale, as in the CMSSM studied here, the (mi/2,mo) planes usually studied for 
fixed tan/3 and sign of /i are effectively reduced to narrow strips of limited thickness in 
mo for any given value of mi/2 [12] and the other parameters. 

These strips have been delineated and parametrized when = for several choices 
of tan/3 for each sign of /x, and the possible LHC and ILC phenomenology along these 
lines has been discussed [16]. As preliminaries to studying indirect sensitivities to the scale 
of supersymmetry along some of these WMAP strips, we first address a couple of physics 
issues. One is that the experimental central value of rut has changed since Ref. [16], from 
174.3 GeV to 178.0 GeV [21], and the other is the dependence of the WMAP strips on Aq. 
The change in rrit has a significant effect on the regions of CMSSM parameter space allowed, 
particularly in the focus-point region where the range of mo allowed by cosmology now starts 
above 4 TeV. In view of the high values of mo and the sensitivity to rrit [22] , we do not study 
the focus-point region further in this paper. There are also rrit- and 74o-dependent effects 
in the 'funnels' where neutrahnos annihilate rapidly via the H, A poles. These affect the 
dependence of mo on mi/2 along the WMAP lines, as we now discuss in more detail. As 
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explained below, because of the anomalous magnetic moment of the muon, we focus on cases 
with /i > 0. 

Plotted in Fig. [T] is the region in the (mi/2,mo) plane for fixed tan/?,Ao and /i > for 
which the relic density is in the WMAP range (the results of [9] are in qualitative agreement 
with Ref. [23]). We have applied cuts based on the lower limit to the Higgs mass, b 57, and 
require that the LSP be a neutralino rather than the stau. The thin strips correspond to the 
relic density being determined by either the coannihilation between nearly degenerate f's and 
x's or, as seen at high tan (3, by rapid annihilation when ^ M4/2. We see in Fig. ^a) that 
the WMAP strip for /i > and tan /5 = 10 does not change much as rrit is varied, reflecting 
the fact that the allowed strip is dominated by annihilation of the neutralino LSP x with the 
lighter stau slepton fi. The main effect of varying rrit is that the truncation at low mi/2, due 
to the Higgs mass constraint, becomes more important at low rrit. This effect is not visible 
in Fig. ^b) for tan/3 = 50, where the cutoff at low mi/2 is due to the b s'j constraint, 
and rapid XX ^ A,H annihilation is important at large mi/2. The allowed regions at larger 
mi/2 vary significantly with rrit when tan/? = 50, because the A,H masses and hence the 
rapid-annihilation regions are very sensitive to rrit through the renormalization group (RG) 
running. Indeed, the rapid-annihilation region almost disappears for rut = 182 GeV at this 
value of tan (3. In this case, in particular, we see a wisp of allowed CMSSM parameter space 
running almost parallel to, but significantly above, the familiar coannihilation strip, which is 
due to rapid fifi H annihilation. At higher values of tan/5 the rapid-annihilation region 
would reappear for rrit = 182 GeV. 

We now turn to the variation of the WMAP strips for different Aq, but with rrit fixed 
to nit = 178 GeV. Since the WMAP strips are largely independent of the sign of /i, for 
clarity we show them in Fig. El only for /i > 0. We see in Fig. |2fa,b) that the WMAP strip 
for tan/5 = 10 also does not change much as Aq is varied: the main effect is for the strip 
to move to larger mo as \Aq\ is increased. This is because the main effect of Aq is on the 
running of the diagonal stau masses, whose RG equations depend only on A^. The splitting 
of the two stau masses depends on the sign of Aq via the off-diagonal entries in the stau 
mass matrix, but the impact of this effect on the final stau masses is relatively small. Hence 
the WMAP strips rise for both signs of Aq. For a given value of mi/2,rnQ and tan/5, the 
low-energy value of Ar is shifted from its high-energy value, Aq, by an amount AA that 
is relatively independent of Aq. Therefore, for \Aq\ much larger than AA, the low-energy 
value of At- will be larger than that for Aq = 0, causing the right-handed stau soft mass to 
drop. This in turn increases the value of mo corresponding to the coannihilation strip. Only 
when the low-energy value of l^,-! is less than and of opposite sign to A A does the light stau 
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Figure 1: The WMAP strips for /i > 0, = and (a) tan/? = 10, (b) tan/5 = 50, showing 
the dependence on the top-quark mass, for rrit = 174, 178 and 182 GeV. 



mass increase. In the specific examples shown in Fig. |2fa,b), AA ranges from about 130 
GeV at low mi/2 to about 550 GeV at high mi/2- Since the shifts are always positive, the 
coannihilation strip rises less for negative values of Aq (Fig. [2fb)) than for positive values 
(Fig. Ha)). 

The WMAP regions for tan/5 = 50 vary much more rapidly with \Aq\, because of the 
sensitivity of the A, H masses and hence the rapid-annihilation regions. In Fig. |21^c) the case 
for v4o > can be seen, whereas Fig. |2l^d) shows < 0. We again see wisps of allowed 
CMSSM parameter space due to rapid fifi H annihilation. In this described 
above, the right-handed stau mass is sensitive to the value of Aq. Therefore, for 7^ 
(Fig. |21^c,d)), the cosmologically preferred region shifts to larger mg for both signs of Aq. In 
addition, the value of the heavy Higgs scalar and pseudoscalar masses depends on Aq (not 
only Aq) and the position of the rapid-annihilation funnels therefore depends sensitively 
on Aq. 

In the following, we mainly present our results along the WMAP strips for = 178 GeV, 
the present experimental central value [21], but we do show results for different values of |Ao|. 
This is because the variation with is less important for tan/5 = 10, and comparable with 
that due to varying |y4o| when tan/5 = 50. Additionally, we present scans of the (mi/2, ^o) 
planes for tan /? = 10 and 50. 
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Figure 2: The WMAP strips for /i > 0, = 178 GeV and (a) tan/? = 10, > (upper 
left), (b) tan /? = 10, < (upper right), (c) tan /? = 50, > (lower left), (d) tan /? = 50, 
< (lower right) showing the dependence on Aq for Aq = 0, ±mi/2 and ±2mi/2- 
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3 Present and future sensitivities to the scale of super- 
symmetry from low-energy observables 



In this section, we briefly describe the low-energy observables used in our analysis. We discuss 
the current and prospective future precision of the experimental results and the theoretical 
predictions. In the following, we refer to the theoretical uncertainties from unknown higher- 
order corrections as 'intrinsic' theoretical uncertainties and to the uncertainties induced by 
the experimental errors of the input parameters as 'parametric' theoretical uncertainties. We 
also give relevant details of the higher-order perturbative corrections that we include. We do 
not discuss theoretical uncertainties from the RG running between the high-scale parameters 
and the weak scale (see Ref. [19] for a recent discussion in the context of predicting the 
CDM density). At present, these uncertainties are expected to be less important than 
the experimental and theoretical uncertainties of the precision observables. In the future, 
both the uncertainties from unknown higher-order terms in the RG running and from the 
parameters entering the running will considerably improve. 

Results for these observables are shown as a function of mi/2 with Aq varied, mo deter- 
mined by the WMAP constraint (see Sect.E)), and tan/5 = 10,50. In this way the indirect 
sensitivities of the low-energy observables to the scale of supersymmetry are investigated. 

3.1 The W boson mass 

The W boson mass can be evaluated from 

where a is the flne structure constant and Gp the Fermi constant. The radiative corrections 
are summarized in the quantity Ar [24]. The prediction for Mw within the Standard Model 
(SM) or the MSSM is obtained from evaluating Ar in these models and solving eq. (0) in an 
iterative way. 

The one-loop contributions to Ar can be written as 

Ar = Aa-4^Ap + (Ar),em, (2) 

where Aa is the shift in the flne structure constant due to the light fermions of the SM, 
Aa oc logruf, and Ap is the leading contribution to the p parameter. It is given by fermion 
and sfermion loop contributions to the transverse parts of the gauge boson self-energies at 
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zero external momentum, 

Sz(0) S„,(0) 

The remainder part, (Ar)rem, contains in particular the contributions from the Higgs sector. 

We include the complete one-loop result in the MSSM [25,26] as well as higher-order QCD 
corrections of SM type of 0{aas) [27,28] and O^aaf) [29,30]. Furthermore, we incorporate 
super symmetric corrections of 0{aas) [31] and of 0{af) [32] to Ap. 

The remaining intrinsic theoretical uncertainty in the prediction for within the 
MSSM is still significantly larger than in the SM, where it is currently estimated to be 
about 4 MeV [33]. We estimate the present [34] and future intrinsic uncertainties to be 

^^intr,current ^ y^^y ^ ^^intr,future ^ ^ y^^y _ 

The parametric uncertainties are dominated by the experimental error of the top-quark mass 
and the hadronic contribution to the shift in the fine structure constant. The current errors 
induce the following parametric uncertainties 

^^current ^ 4 3 g^y ^ ^^para,m„current _ ^ 

5(Aa™*) = 36 X 10"^ aMP?'""'^"''^'^'™"^"* ^ 6.5 MeV. (6) 

At the ILC, the top-quark mass will be measured with an accuracy of about 100 MeV [2]. 
The parametric uncertainties induced by the future experimental errors of rrit and Aahad [35] 
will then be [36] 

5{Aa^^^'') = 5 X 10-^ AMP,"=""'^°'-<^'f"*"^ ^ 1 MeV. (8) 
The present experimental value of is [1] 

^^P,CUrrent _ ^^^^^^ ^ q gg^ ^^^y^ ^g^ 

With the GigaZ option of the ILC (i.e. high-luminosity running at the Z resonance and the 
WW threshold) the VT-boson mass will be determined with an accuracy of about [37, 38] 

In all plots of this section we show the theory predictions without parametric and intrinsic 
theoretical uncertainties (using nit = 178 GeV). In the fits carried out in Sects. IHandElbelow 
we take both parametric and intrinsic theoretical uncertainties into account. 
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Figure 3: The CMSSM prediction for M]y as a function of mi/2 along the WMAP strips for 
(a) tan/? = 10 and (b) tan/3 = 50 for various Aq values. In each panel, the centre (solid) 
line is the present central experimental value, and the (solid) outer lines show the current 
±1-0" range. The dashed lines correspond to the anticipated GigaZ accuracy, assuming the 
same central value. 



We display in Fig. IHl the CMSSM prediction for Mw and compare it with the present 
measurement (solid lines) and a possible future determination with GigaZ (dashed lines). 
Panel (a) shows the values of obtained with tan/? = 10 and |y4o| < 2, and panel (b) 
shows the same for tan/5 = 50. It is striking that the present central value of Mvi/ (for 
both values of tan/?) favours low values of mi/2 ~ 200-300 GeV, though values as large as 
800 GeV are allowed at the 1-a level, and essentially all values of mi/2 are allowed at the 
90% confidence level. The GigaZ determination of Mw might be able to determine indirectly 
a low value of mi/2 with an accuracy of ±50 GeV, but even the GigaZ precision would still 
be insufficient to determine mi/2 accurately if mi/2 ^ 600 GeV. 

3.2 The effective leptonic weak mixing angle 

The effective leptonic weak mixing angle at the Z boson resonance can be written as 

sin'e^, = ] (l-Re^) , (11) 

where v^s and a^s denote the effective vector and axial couplings of the Z boson to charged 
leptons. As in the case of Mw, the leading supersymmetric higher-order corrections enter 
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via the p parameter, 

5sin2^eff^-#^Ap. (12) 

Our theoretical prediction for sin^ contains the same higher-order corrections as described 
in Sect. 13.11 

In the SM, the remaining intrinsic theoretical uncertainty in the prediction for sin^ ^eff 
has been estimated to be about 5 x 10~^ [39]. For the MSSM, we use as present [34] and 
future intrinsic uncertainties 

A sin^ = 12 X 10-^ A sin^ ^'I'^.f"*-- = 2 x 10"^ (13) 

The current experimental errors of and Aa^ad induce the following parametric uncertain- 
ties 

5m™* = 4.3 GeV ^ A sin^ ^P|-.-t.™t ~ 14 x 10-^ (14) 

5(Aa™*) = 36 X 10-^ ^ ^ gi^2 ^para,Aahad,currcnt _ ^3 ^ ^q-5_ ^^5) 

These should improve in the future to 

^^futurc = 0.1 GeV ^ A Sin^ ^P--.-*.f"t-e _ ^ ^q-S^ ^^g) 

5(A4"X') = 5 X 10-^ ^ Asin^^P^'^^'^^'^-'^'f"*"^ ^ 1.8 X 10"^ (17) 

It is well known that there is a 2.8-cr discrepancy [1] between the leptonic and heavy- 
flavour determinations of the electroweak mixing angle, with the leptonic measurement of 
sin^ 6'efr tending to pull down the value of Higgs-boson mass preferred in the SM fit, whereas 
the heavy-flavour measurements favour a larger value of the Higgs mass. The Electroweak 
Working Group notes that the overall quality of a global electroweak flt is quite acceptable, 
~ 26% [1], and we use their combination of the two sets of measurements: 

sin' ^exp.current ^ 0.23150 ± 0.00016. (18) 

The experimental accuracy will improve to about 

^sin'^^^P'f"*"--^ = 1 X 10-^ (19) 

at GigaZ [40]. 

Fig. m shows the prediction for sin' in the CMSSM compared with the present and 
future experimental precision. As in the case of Mw-, low values of mi/2 are also favoured 



10 



0.2320 



0.2315 



0.2310 



0.2305 



0.2300 



CMSSM, \i>0 
tanp = 10, A„ = 
tanp = 10, A„ = +m^,2 
tanp = 10, A„ = -m,,2 
tanp = 10, A^ = +2 
tanp = 10, A = -2 m 



200 



400 



600 



800 



1000 1200 1400 



0.2320 



0.2315 



0.2310 



0.2305 



0.2300 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


«. -8- -•c* ? 




----B--^ 

!S8S° 






.•: 


CMSSM, |j > 
. tanp = 50, A„ = 




• tanp = 50, A„ = +m^,2 
tanp = 50, A„ = -m^,^ 
. tanp = 50, A„ = +2 m,,^ 


, < 1 , 


o tanp = 50, A„ = -2 m,,^ 
, 1 < , , 1 , , < 1 , < , 1 < , , 1 , , < 1 , 



200 



400 



600 



800 1000 1200 1400 



m,,, [GeV] 



m, [GeV] 



Figure 4: The CMSSM prediction for sin^^cfj as a function of mi/2 along the WMAP strips 
for (a) tan/3 = 10 and (b) tan/5 = 50 for various Aq values. In each panel, the centre (solid) 
line is the present central experimental value, and the (solid) outer lines show the current 
±l-cr range. The dashed lines correspond to the anticipated GigaZ accuracy, assuming the 
same central value. 



independently by sin^^efr- The present central value prefers mi/2 = 300-500 GeV, but the 
1-cr range extends beyond 1500 GeV (depending on Aq), and all values of mi/2 are allowed at 
the 90% confidence level. The GigaZ precision on sin^ 6cs would be able to determine mi/2 
indirectly with even greater accuracy than Mw at low mi/2, but would also be insufficient 
if mi/2 > 700 GeV. 

3.3 The anomalous magnetic moment of the muon 

We now discuss the evaluation of the MSSM contributions to the anomalous magnetic mo- 
ment of the muon, = {g — 2)^. Since the possible deviation of the SM prediction from 
the experimental result is crucial for the interpretation of the results, we first review this 
aspect in the light of recent developments. 

The SM prediction for the anomalous magnetic moment of the muon (see Refs. [41,42] for 
reviews) depends on the evaluation of the hadronic vacuum polarization and light-by-light 
(LBL) contributions. The former have been evaluated in [43-46] and the latter in [47,48]. 
The evaluations of the hadronic vacuum polarization contributions using e~^e~ and r decay 
data give somewhat different results. Recently, new data have been published by the KLOE 
Collaboration [49], which agree well with the previous data from CMD-2. This, coupled 
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with a greater respect for the uncertainties inherent in the isospin transformation from r 
decay, has led to a proposal to use the e"'"e~ alone and shelve the r data, resulting in the 
estimate [50] 

af''° = (11 659 182.8 ± 6.3had ± 3.5lbl ± 0.3qed+ew) x 10"!°, (20) 

where the source of each error is labelled ^. 

This result is to be compared with the final result of the Brookhaven {g 
E821, namely [52] 

a;i^P = (11 659 208.0 ± 5.8) x 10"^°, 
leading to an estimated discrepancy 

a^"P - = (25.2 ± 9.2) x 10"^°, (22) 

equivalent to a 2.7 a effect. In view of the chequered history of the SM prediction, eq. (j20|l . 
and the residual questions concerning the use of the r decay data, it would be premature 
to regard this discrepancy as firm evidence of new physics. We do note, on the other hand, 
that the ((? — 2)^ measurement imposes an important constraint on supersymmetry, even if 
one uses the r decay data. We use eq. for our numerical discussion below. 

The following MSSM contributions to the theoretical prediction for have been con- 
sidered. We take fully into account the complete one- loop contribution to a^, which was 
evaluated nearly a decade ago in Ref. [53]. We make no simplification in the sparticle mass 
scales but, for illustrating the possible size of corrections, a simplified formula can be used, 
in which relevant supersymmetric mass scales are set to a common value, MgusY = rny.± = 
m^o = rrifi = nii,^. The result in this approximation is given by 

susY.iL ^ 13 ^ iQ-10 (^M^f tan/? sign(/i). (23) 
^ V MsuSY ) ^ ^ > 

We see that supersymmetric effects can easily account for a (20 . . . 30) x 10^^° deviation, if /x 
is positive and MgusY lies roughly between 100 GeV (for small tan/5) and 600 GeV (for large 
tan (3). For this reason, in the rest of this paper, we restrict our attention to ;U > 0. Even in 
view of the possible size of experimental and theoretical uncertainties, it is very difficult to 
reconcile /i < with the present data on a^. 

In addition to the full one-loop contributions, we also include several two-loop corrections. 
The first class of corrections comprises the leading log (^^/Msusy) terms of supersymmetric 
"'^The updated QED result from [51] is included. 
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one-loop diagrams with a photon in the second loop, which are given by [54]: 



(24) 



These amount to about —8% of the supersymmetric one-loop contribution for a supersym- 
metric mass scale Msusy = 500 GeV. 

The second class of two-loop corrections comprises diagrams with a closed loop of SM 
fermions or scalar fermions. These were calculated in Ref. [55], where it was demonstrated 
that these corrections may amount to ~ 5 x 10~^° in the general MSSM, if all experi- 
mental bounds are taken into account. These corrections are included in the Fortran code 
FeynHiggs [56,57]. We have furthermore taken into account the 2-loop contributions to 

from diagrams containing a closed chargino/neutralino loop, which have been evalu- 
ated in [58]. Here we use an approximate form for these corrections, which are typically 
~ 1 X 10"i°. 

The current intrinsic uncertainties in the MSSM contributions to can be estimated 
to be ~ 6 X 10~^° [58,59]. In the more restricted CMSSM parameter space the intrinsic 
uncertainties are smaller, being about 1 X 10-1°. Once the full two-loop result in the MSSM 
is available, this uncertainty will be further reduced. We assume that in the future the 
uncertainty in eq. will be reduced by a factor two. 

As seen in Fig. the CMSSM prediction for is almost independent of Aq for tan j3 = 
10, but substantial variations are possible for tan/? = 50, except at very large mi/2- In the 
case tan/3 = 10, mi/2 ~ 200-400 GeV is again favoured at the ±1-0" level, but this preferred 
range shifts up to 400 to 800 GeV if tan/? = 50, depending on the value of Aq. At the 2-a 
level, there is nominally an upper bound mi/2 ^ 600(1100) GeV for tan/3 = 10(50), but 
according to the discussion above it should be interpreted with care. Nevertheless, the lower 
bound to mi/2 for both tan/? = 10 and 50 should be regarded as relatively robust. On the 
other hand, it is striking that M^y, sin^ ^efr and all favour small mi/2 for tan/3 = 10. If 
tan /3 = 50, the consistency between the ranges preferred by the different observables is not 
so striking. 

3.4 The decay 6 ^ 57 

Since this decay occurs at the loop level in the SM, the MSSM contribution might, a priori, 
be of similar magnitude. The most up-to-date theoretical estimate of the SM contribution 
to the branching ratio is [60] 



BR(6 ^ 57) = (3.70 ± 0.30) x 10 



,-4 



(25) 



13 



60 



50 



40 - 



o 

ZZ. 30 

CO 

20 



10 



o • 
• o « 



CMSSM,^>0 

• tanp = 10 

• tanp = 50, \ = 

• tanp = 50, Aq = 
tanp = 50, Ap = -m^^ 



m tanp = 50, Aq = +2 
o tanp = 50, Aq = -2 



/2 




200 



400 



600 



800 1000 1200 1400 



[GeV] 



Figure 5: The CMSSM prediction for Aa^ as a function of mi/2 along the WMAP strips for 
tan (3 = 10, 50 and different Aq values. The central (solid) line is the central value of the 
present discrepancy between experiment and the SM value evaluated using e~^e~ data (see 
text), and the other solid (dotted) lines show the current ±l(2)-cr ranges, eq. 



where the calculations have been carried out completely to NLO in the MS renormalization 
scheme, and the error is dominated by higher-order QCD uncertainties. A complete NNLO 
QCD calculation is now underway, and will reduce significantly the uncertainty, once it is 
available. 

For comparison, the present experimental value estimated by the Heavy Flavour Aver- 
aging Group (HFAG) is [61] 

BR(6 ^ S7) = (3.54^0;^) x 10-^ (26) 

where the error includes an uncertainty due to the decay spectrum, as well as the statistical 
error. The very good agreement between eq. (j26|) and the SM calculation eq. (|25|) imposes 
important constraints on the MSSM, as we see below. 

Our numerical results have been derived and checked with three different codes. The 
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first is based on Refs. [62,63]^ and the second is based on Refs. [63,64]^. Results have been 
derived using the charm pole mass as well as the charm running mass, giving an estimate 
of remaining higher-order uncertainties. Finally, our results have been checked with the 
BR(6 — s> 57) evaluation provided in Ref. [65], which yielded very similar results to our 
two other approaches. For the current theoretical uncertainty of the MSSM prediction for 
BR(6 — » 57) we use the value of eq. ()25|) . For the future uncertainty from the experimental 
as well as the theoretical side we assume a reduction by a factor of 3. 
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Figure 6: The CMSSM predictions for BR(6 — >■ 57) as a function of mi/2 along the WMAP 
strips for (a) tan /? = 10 and (b) tan /3 = 50 and various choices of Aq. The uncertainty shown 
combines linearly the current experimental error and the present theoretical uncertainty in 
the SM prediction. The central (solid) line indicates the current experimental central value, 
and the other solid (dotted) lines show the current ±1(2)-(T ranges. The dash-dotted line 
corresponds to a more conservative estimate of intrinsic uncertainties (see text). 



As already mentioned, the present central value of this branching ratio agrees very well 

with the SM, implying that large values of mi/2 cannot be excluded for any value of tan/5. 

The uncertainty range shown in Fig. EJ combines linearly the current experimental error 

and the present theoretical uncertainty in the SM prediction. Note however, that at present 

there is also an uncertainty in the computed MSSM value (included in obtaining the excluded 

regions in Figs. Q and ^ from the uncertainty in the SUSY loop calculations. Taking this 

conservatively into account results in a 95% CL. exclusion bound of 0.00016 in the case of 

tan P = 10, and of 0.000195 in the case of tan /3 = 50. These values are shown as dash-dotted 

^We are grateful to P. Gambino and G. Ganis for providing the corresponding code. 
■^We thank Gudrun HiUer for providing the corresponding Fortran code. 
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lines in Fig. IHl This allows a somewhat lower range in mi/2 than depicted in Fig. IHl We 
assume that these uncertainties can be significantly reduced in the future. We have checked 
that they have no significant impact on the results presented below. 

Since the CMSSM corrections are generally smaller for smaller tan/3, even values of mi/2 
as low as ~ 200 GeV would be allowed at the 90% confidence level if tan/5 = 10, whereas 
~ 450 GeV would be required if tan/? = 50. These limits are very sensitive to Aq, and, 
if the future error in BR(6 — > 57) could indeed be reduced by a factor ~ 3, the combination 
of BR(6 — > 57) with the other precision observables might be able, in principle, to constrain 
Aq significantly. 

3.5 The branching ratio Bg ii^fi~ 

The SM prediction for this branching ratio is (3.4 ± 0.5) x 10^^ [66], and the present exper- 
imental upper limit from the Fermilab Tevatron collider is 3.4 x 10"'' at the 95% CL. [67], 
providing ample room for the MSSM to dominate the SM contribution. The current Tevatron 
sensitivity, being based on an integrated luminosity of about 410 pb^^ summed over both 
detectors, is expected to improve significantly in the future. A naive scaling of the present 
bound with the square root of the luminosity yields a sensitivity at the end of Run II of 
about 5.4 X 10~^ assuming 8 fb~^ collected with each detector. An even bigger improvement 
may be possible with better signal acceptance and more efficient background reduction. In 
Ref. [68] an estimate of the future Tevatron sensitivity of 2 x 10~^ at the 90% CL. has been 
given, and a sensitivity even down to the SM value can be expected at the LHC Assuming 
the SM value, i.e. BR(i?s — > f^^f^ ) ~ 3.4 x 10~^, it has been estimated [69] that LHCb can 
observe 33 signal events over 10 background events within 3 years of low-luminosity running. 
Therefore this process offers good prospects for probing the MSSM. 

For the theoretical prediction we use results from Ref. [70]^, which include the full one- 
loop evaluation and the leading two-loop QCD corrections. We are not aware of a detailed 
estimate of the theoretical uncertainties from unknown higher-order corrections. 

In Fig. Elthe CMSSM prediction for BR(i?s —>■ fi'^fj^~) as a function of mi/2 is compared 
with the present Tevatron limit and our estimate for the sensitivity at the end of Run II. 
For tan /3 = 10 the CMSSM prediction is significantly below the present and future Tevatron 
sensitivity. With the current sensitivity, the Tevatron starts to probe the CMSSM region 
with tan /? = 50. The sensitivity at the end of Run II will test the CMSSM parameter space 
with tan/3 = 50 and mi/2 ^ 600 GeV, in particular for positive values of Aq. The LHC will 
''We are grateful to A. Dedes for providing the corresponding code. 
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Figure 7: The CMSSM prediction for BR(i?s as a function of mi/2 along the 

WMAP strips for tan /3 = 10 and all Aq values, and for tan /5 = 50 with various values of 
Aq. The solid line shows the current Tevatron limit at the 95% C.L., and the dotted line 
corresponds to an estimate for the sensitivity of the Tevatron at the end of Run II. 



be able to probe the whole CMSSM parameter space via this rare decay. 



3.6 The lightest MSSM Higgs boson mass 

The mass of the lightest CP-even MSSM Higgs boson can be predicted in terms of the other 
CMSSM parameters. At the tree level, the two CP-even Higgs boson masses are obtained as a 
function of M^, the CP-odd Higgs boson mass Ma-, and tan (3. In the Feynman-diagrammatic 
(FD) approach, which we employ here, the higher-order corrected Higgs boson masses are 
derived by finding the poles of the iJ-propagator matrix. This is equivalent to solving 



P 



m 



/i,tree 



^hh{p^ 



X 



P 



m 



H,tT:ee 



^hH{p^ 







(27) 



where the S(p^) denote the renormahzed Higgs-boson self-energies, and p is the external 
momentum. 
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For the theoretical prediction of Mh we use the code FeynHiggs [56,57], which includes 
all numerically relevant known higher-order corrections. The status of the incorporated 
results for the self-energy contributions to eq. ()27|) can be summarized as follows. For the 
one-loop part, the complete result within the MSSM is known [71-73]. Concerning the 
two-loop effects, their computation is quite advanced, see Ref. [74] and references therein. 
They include the strong corrections at (9(at«s) and Yukawa corrections at 0{af), as well 
as the dominant one-loop 0{at) term, and the strong corrections from the bottom/sbottom 
sector at 0{ab<ys)- For the b/b sector corrections also an all-order resummation of the tan jS - 
enhanced terms, tan/5)"), is known [75,76]. Most recently, the 0{atab) and 0{al) 

corrections have been derived [77]. ^ 

The current intrinsic error of due to unknown higher-order corrections and its prospec- 
tive improvement in the future have been estimated to be [74, 79] 

^^mtr,current ^ 3 ^ ^^intr,futurc _ q 5 ^28) 

The estimated future uncertainty assumes that a full two-loop result, leading three- loop and 
possibly even higher-order corrections become available. 

Concerning the parametric error on M^, the top-quark mass has the largest impact, 
entering oc mf at the one-loop level. As a rule of thumb, an uncertainty of drrit = 1 GeV 
translates to an induced parametric uncertainty in Mh of AM™' ^ 1 GeV [80]. We find for 
the parametric uncertainties induced by the present experimental errors of nit and a.^ 

5m™* = 4.3 GeV ^^para,™,, current _ ^ ^ ^29) 

= 0.002 ^ ^J^^para,a„current _ q 3 ^^^^ 

These will improve in the future to 

^^futurc ^ Q^y ^ ^^para.m,, future _ ^^y ^ ^3^^ 

= 0.001 AMr"'"''^"'"" ~ 0.1 GeV. (32) 

Thus, the intrinsic error would be the dominant source of uncertainty in the future. On the 
other hand, a further reduction of the unknown higher-order corrections to Mh is in principle 
possible. 

The experimental accuracy on M^ at the ILC [2] will be even higher than the prospective 
precision of the theory prediction, 

^^exp,future _ q Q^y ^33^ 



^Furthermore, a two-loop effective potential calculation has been carried out in Ref. [78], but no public 
code based on this result is available. 



18 



CMSSM, n > 
tanp= 10, Aj, = 
tanp= 10, Aj, = +m^,2 
tanp =10, A(, = -m,,2 
tanp= 10, Aj, = +2m,^ 
tanp = 10, A„ = -2 m,,. 



130 



120 



> 
C5 



200 



400 



600 800 
%2\ 



1000 1200 1400 



100 





.- • 

- ,m-- 


1 1 1 1 1 

oooooooooa 


• • 

• • 

• 


'„ ocP°°° o ° 
o o o° 
oo 






* 






• . V °° 






* . o" 
•.Vo°o° 








CMSSM, > 
tanp = 50, Aj, = 





• 


tanp = 50, A„ = 


+m„2 


o 


tanp = 50, A„ = 




• 


tanp = 50, A^ = 


+2 


o 


tanp = 50, A„ = 

1 , 1 1 , , 


-2 



200 



400 



600 



800 



1000 1200 



1400 



m,„ [GeV] 



Figure 8: The CMSSM predictions for Mh as functions of mi/2 with (a) tan/? = 10 and 
(b) tan/? = 50 for various Aq. A hypothetical experimental value is shown, namely Mh = 
120 GeV. We display an optimistic anticipated theory uncertainty of ±0.2 GeV, as well 
as a more realistic theory uncertainty of ±0.5 GeV and the current theory uncertainty of 
±3 GeV. 



We show in Fig. |Hlwe show the for M^, assuming a hypothetical measurement at Mh = 
120 GeV. Since the experimental error at the ILC will be smaller than the prospective 
theory uncertainties, we display the effect of the current and future intrinsic uncertainties. In 
addition, a more optimistic value of 200 MeV is also shown. The figure clearly illustrates the 
high sensitivity of this electroweak precision observable to variations of the supersymmetric 
parameters (detailed results for Higgs boson phenomenology in the CMSSM can be found in 
Ref. [81]). The comparison between the measured value of Mh and a precise theory prediction 
will allow one to set tight constraints on the allowed parameter space of mi/2 and Aq. 

3.7 The Higgs boson branching ratios 

Within the CMSSM, various Higgs boson decay channels will be accessible at the LHC and 
the ILC. At the LHC, Higgs boson couplings [82] or ratios of them [83,84] can in general be 
determined at the level of ~ 10% at best, depending on the Higgs-boson mass and theoretical 
assumptions. Therefore we concentrate on ILC measurements and accuracies. 
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It has been shown in Ref. [85] that the observable combination 



BR{h - 


bb)/BR{h WW*) 


MSSM 


BR{h - 


bb)/BR{h WW*) 


SM 



(34) 



of Higgs boson decay rates is particularly sensitive to deviations of the MSSM Higgs sector 
from the SM. Even though the experimental error on the ratio of the two branching ratios 
is larger than that on the individual ones, the quantity r has a stronger sensitivity to Ma 
than any single branching ratio. 

For the evaluation of BR{h bb), we use the results of Ref. [86], including the result of 
resumming the contributions of 0((as tan/?)") [75,76]. The evaluation of BR(/i WW*) is 
based on an effective-coupling approach, taking into account off-shell effects. The corrections 
used for the effective-coupling calculation are the same as for the Higgs-boson mass calcula- 
tion, including the full one-loop and leading and subleading two-loop contributions [56,74]. 
The evaluation has been performed with FeynHiggs [56,57]. 

For the prospective accuracy at the ILC, we consider two cases. At the ILC with ^/s = 
500 GeV an accuracy of 4% seems to be feasible [2], whilst at y/s = 1 TeV this accuracy 
could be improved to [87] 

Cr \ exp, future 
= 1.5%. (35) 

Since in this ratio of branching ratios many theoretical uncertainties cancel, we assume that 
the future theoretical error can be neglected. In the analysis in Sect. El we use the accuracy 
of eq. 

In Fig. Elthe results for r are shown as functions of mi/2 for tan/5 = 10, 50. In the figure 
we indicate accuracies of both 4% and 1.5%. For low tan/5, the high ILC accuracy in r will 
allow one to detect a deviation from the SM prediction for all CMSSM points. For large 
tan (3, the effects of the supersymmetric contributions to r are in general smaller. Deviations 
up to mi/2 ~ 1 TeV could be visible, depending somewhat on Aq. 

4 Combined Sensitivity: Present Situation 
4.1 Best fits for WMAP strips at fixed 

We now investigate the combined sensitivity of the four low-energy observables for which 
experimental measurements exist at present, namely Myy, sin^^efr, {9~'^)n ^^^d BR{b sj). 
Since only an upper bound exists for BR(i?s fi^fi^), we discuss it separately below. We 
begin with an analysis of the sensitivity to mi/2 moving along the WMAP strips with fixed 
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Figure 9: The CMSSM predictions for [BR{h bb)/BR{h -> WW*)]MssM/[BR{h -> 
bb)/BR{h WW*)]sM as functions of mi/2 for (a) tan/? = 10 and (b) tan/3 = 50 with 
various values of Aq. The central (solid) hne corresponds to the SM expectation. The outer 
(dotted) and inner (solid) lines indicate an ILC measurement with 4% and 1.5% accuracy, 
respectively. 



values of Ao and tan/5. The experimental central values, the present experimental errors 
and theoretical uncertainties are as described in Sect. 01 The experimental uncertainties, 
the intrinsic errors from unknown higher-order corrections and the parametric uncertainties 
have been added quadratically, except for BR(6 57), where they have been added linearly. 
Assuming that the four observables are uncorrelated, a fit has been performed with 

A'^ / pexp ptheo \ ^ 

X'^T.(- ^ • (36) 

n=l V / 

Here denotes the experimental central value of the nth observable, so that = 4 for 
the set of observables included in this fit, is the corresponding CMSSM prediction 

and (T„ denotes the combined error, as specified above. We have rejected all points of the 
CMSSM parameter space with either < 113 GeV [88,89] or a chargino mass lighter than 
103 GeV [90]. 

The results are shown in Fig. 1101 for tan/5 = 10 and tan/? = 50. They indicate that, 
already at the present level of experimental accuracies, the electroweak precision observables 
combined with the WMAP constraint provide a sensitive probe of the CMSSM, yielding 
interesting information about its parameter space. For tan/5 = 10, the CMSSM provides a 
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Figure 10: The results of fits based on the current experimental results for the precision 
observables Mw, sin^^efr, {g — 2)^ and BR(6 — 57) are shown as functions of mi/2 in the 
CMSSM parameter space with CDM constraints for different values of Aq. The upper plot 
shows the results for tan /3 = 10, and the lower plot shows the case tan j3 = 50. 
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very good description of the data, resulting in a remarkably small minimum value. The 
fit shows a clear preference for relatively small values of mi/2, with a best-fit value of about 
mi/2 = 300 GeV. The best fit is obtained for Aq < 0, while positive values of Ao result in a 
somewhat lower fit quality. The fit yields an upper bound on mi/2 of about 600 GeV at the 
90% C.L. (corresponding to Ax^ < 4.61). 

These results can easily be understood from the analysis in Sect. 01 For tan/5 = 10, the 
CMSSM prediction with mi/2 ~ 300 GeV is very close to the experimental central values 
of M]Y, sin^6'eff and {g — 2)^ for all values of Aq, see Figs. OHSl Also, BR(6 — > 57) is well 
described for mi/2 ~ 300 GeV and Aq < 0, while large positive values of Aq lead to a CMSSM 
prediction for BR(6 — * 57) which is significantly below the experimental value. Consequently, 
in the case of tan /3 = 10, a very good fit quality is obtained for mi/2 ~ 300 GeV and Aq < 0.^ 
Some of the principal contributions to the increase in when mi/2 increases for tan/5 = 10 
are as follows. For Aq = — ^^1/2, ""^1/2 = 900 GeV, we find that {g — 2)^ contributes about 
5 to Ax^, M\Y nearly 1 and s\v? 0^^ about 0.2, whereas the contribution of BR(6 S7) is 
negligible. On the other hand, for Aq = +2mi/2, which is disfavoured for tan/5 = 10, the 
minimum in is due to a combination of the four observables, but {g — 2)^ again gives the 
largest contribution for large mi/2. 

For tan/5 = 50 the overall fit quality is worse than for tan/5 = 10, and the sensitivity to 
mi/2 from the precision observables is lower. This is related to the fact that, whereas 
and sin^^efj prefer small values of mi/2 also for tan/5 = 50, as seen in Figs. El and EJ the 
CMSSM predictions for {g — 2)^ and BR(6 — > S7) for high tan j3 are in better agreement 
with the data for larger mi/2 values, as seen in Figs. El and El Also in this case the best fit is 
obtained for negative values of Aq, but the preferred values for mi/2 are 200-300 GeV higher 
than for tan/5 = 10. 

In Figs. ITTVIT^ the fit results of Fig. ^1 are expressed in terms of the masses of different 

supersymmetric particles. Fig. ^2 shows that for tan (3 = 10 the best fit is obtained if the 

lightest supersymmetric particle (LSP), which within the CMSSM is the lightest neutralino, 

is lighter than about 200 GeV (with a best-fit value ~ 100 GeV). The best-fit values for the 

masses of the lighter chargino, the second-lightest neutralino (recall also that m-+ ^ '^x"-'' 

both sleptons and the lighter stau are all below 250 GeV, while the preferred region of the 

masses of the heavier chargino and the heavier neutralinos is about 400 GeV. These masses 

offer good prospects of direct sparticle detection at both the ILC and the LHC. There are 

also some prospects for detecting the associated production of charginos and neutralinos at 

preference for relatively small values of mi/2 within the CMSSM has also been noticed in Ref. [7], 
where only (5 — 2)^ and BR(6 — s- 57) had been analyzed. 
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Figure 11: The contours in the CMSSM with tan/3 = 10 for different sparticle masses, 
based on the fits to the parameter space shown in Fig. ^| The first row shows (left) the 
mass of the neutrahno LSP, m^o, and (right) the mass of the hghter chargino, 
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The second row shows (left) the mass of the heavier chargino, k, m^o, and (right) the 
mass of the lighter stau, nif^ . The selectron masses are shown in tfie third row. 
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Figure 12: The contours in the CMSSM with tan/5 = 50 for different sparticle masses, 
based on the fits to the parameter space shown in Fig. ^| The first row shows (left) the mass 
of the hghtest neutrahno, m^o, and (right) the mass of the hghter chargino, ~ ''^x"' 
The second row shows (left) the mass of the heavier chargino, ?a '^xs' (right) the 
mass of the lighter stau, . The selectron masses are shown in tfie third row. 
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Figure 13: The contours in the CMSSM with tan/5 = 10 for different sparticle masses, 
based on the fits to the parameter space shown in Fig. ^1 The first row shows the scalar 
top masses, m^-^, m^-^. The second row shows the scalar bottom masses, m^,^, m^^. The third 
row shows the gluino mass, m^, (left) and the mass of the scalar Higgs boson. Ma (right). 
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Figure 14: The contours in the CMSSM with tan/3 = 50 for different sparticle masses, 
based on the fits to the parameter space shown in Fig. ^1 The first row shows the scalar 
top masses, m^-^, m^-^. The second row shows the scalar bottom masses, m^,^, m^^. The third 
row shows the gluino mass, m^, (left) and the mass of the scalar Higgs boson. Ma (right). 
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the Tevatron collider, via their trilepton decay signature, in particular. This is estimated to 
be sensitive to mi/2 ~ 250 GeV [91], covering much of the region below the best-fit value of 
mi/2 that we find for tan/5 = 10. 

The same particle masses in the case tan /? = 50 are shown in Fig. ^| Here the best-fit 
values for the LSP mass and the lighter stau are still below about 250 GeV. The minimum 

for the other masses is shifted upwards compared to the case with tan/3 = 10. The best- 
fit values are obtained in the region 400-600 GeV. Correspondingly, these sparticles would 
be harder to detect. At the ILC with ^/s ^ 1 TeV, the best prospects would be for the 
production of X1X2 of ^1^1 • Other particles can only be produced if they turn out to be 
on the light side of the function. 

In Fig. El we focus on the coloured part of the supersymmetric spectrum and the 
Higgs mass scale. The case of tan /? = 10 is shown in Fig. The top row shows the two 
scalar top masses, the middle row displays the two scalar bottom masses, and the bottom 
row depicts the gluino mass and Ma- All the coloured particles should be accessible at the 
LHC. However, among them, only ti has a substantial part of its x^-favoured spectrum below 
500 GeV, which would allow its detection at the ILC. The same applies for the mass of the 
A boson. The Tevatron collider has a sensitivity to m^^ ^ 450GeV, which is not far below 
our best-fit value for tan/? = 10 [91]. 

Finally, in Fig. El we show the same masses in the case of tan/5 = 50. All the particles 
are mostly inaccessible at the ILC, though the LHC has good prospects. However, at the 
90% C.L. the coloured sparticle masses might even exceed ~ 3 TeV, which would render 
their detection difficult. Concerning the heavy Higgs bosons, their masses may well be 
below ~ 1 TeV. In the case of large tan /?, this might allow their detection via the process 
bb blH/A bb t+t- [92]. 

4.2 Scan of the CMSSM Parameter Space 

Whereas in the previous section we presented fits keeping Ao/mi/2 fixed, we now analyse the 
combined sensitivity of the precision observables Mw, sin^ 6'efr, BR(6 — > 57) and {g — 2)^ in a 
scan over the (mi/2,^0) parameter plane. In order to perform this scan, we have evaluated 
the observables for a finite grid in the (^^1/2, "^o) parameter space, fixing mo using the 
WMAP constraint. As before, we have considered the two cases tan/5 = 10 and tan/5 = 50. 
Due to the finite grid size, very thin lines in the (mi/2, ^0) plane for tan/5 = 50, see Fig. |21 
can either be missed completely, or may be represented by only a few points. 

Fig. CH shows the WMAP-allowed regions in the (mi/2,^0) plane for tan/5 = 10 and 
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Figure 15: The results of fits for tan P = 10 (upper plot) and tan P = 50 (lower plot) based 
on the current experimental results for the precision observables Mw, sin^ 6'cfT, {g — 2)^ and 
BR(6 — >■ 57) are shown in the (^1/2, ^0) planes of the CMSSM with the WMAP constraint. 
The best-fit points are indicated, and the coloured regions correspond to the 68% and 90% 
CL. regions, respectively. 
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tan [5 = 50. The current best-fit values obtained via fits for tan (3 = 10 and tan (3 = 50 
are indicated. Tfie coloured regions around tlie best-fit values correspond to the 68% and 
90% C.L. regions (corresponding to Ax^ < 2.30,4.61, respectively). 

For tan/5 = 10 (upper plot of Fig. [T3j). the precision data yield sensitive constraints on 
the available parameter space for mi/2 within the WMAP-allowed region. The precision data 
are less sensitive to Aq. The 90% C.L. region contains all the WMAP-allowed Aq values in 
this region of mi/2 values. As expected from the discussion above, the best fit is obtained 
for negative Aq and relatively small values of mi/2. At the 68% C.L., the fit yields an upper 
bound on mi/2 of about 450 GeV. This bound is weakened to about 600 GeV at the 90% 
C.L. 

As discussed above, the overall fit quality is worse for tan/5 = 50, and the sensitivity 
to mi/2 is less pronounced. This is demonstrated in the lower plot of Fig. ^[ which shows 
the result of the fit in the (mi/2,^0) plane for tan/5 = 50. The best fit is obtained for 
mi/2 ~ 500 GeV and negative Aq. The upper bound on mi/2 increases to nearly 1 TeV at 
the 68% C.L. 

The holes in the coverage of the (mi/2, ^0) plane arise from the finite grid size of the 
scanning procedure, as mentioned above. They would be filled if our scan would also pick up 
the very thin lines, especially the wisps arising from fifi H. Thus, the holes correspond 
to an extremely fine-tuned part of the parameter space, and are sparsely populated but not 
empty. 

In Fig. El we analyze the prospects for the Tevatron to observe the process Bs —>■ /i'^/U". 
We show the regions of the parameter space that are favoured at the 68% or 90% C.L., as a 
result of our fits to the precision observables described above for tan /5 = 10 and tan /? = 50. 
The dotted line corresponds to our estimate of the final Tevatron sensitivity at the 95% C.L. 
of 5.4 X 10~^, see Sect. 13.51 It can be seen that, even for tan/5 = 50, all parameter points 
result in a prediction for BR(i?s /W^/U ) that is below our estimate of the future Tevatron 
sensitivity at the 95% C.L. Only with the more optimistic estimate of 2 x 10^® at the 90% 
C.L., discussed above, could a part of the favoured region for tan/5 = 50 be probed. The 
LHC, on the other hand, will cover the whole CMSSM parameter space. 

5 Combined Sensitivity: ILC Precision 

5.1 Best fits for WMAP strips at fixed ^0 

We now turn to the analysis of the future sensitivities of the precision observables, based on 
the prospective experimental accuracies at the ILC and the estimates of future theoretical 
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Figure 16: Predictions for BR(i?s /^^A* ) within the CMSSM with WMAP constraints are 
shown as functions of mi/2, corresponding to the best-fit regions obtained by a fit (see 
Fig. ITH|l based on the current experimental results for the precision observables My/-, sin^ 6'efr, 
{g — 2)^ and BR(6 —> 37). The different colours indicate the 68% and 90% CL. regions. The 
present bound on BR(i?s /i"*"/!") from the Tevatron (solid line) and our estimate for the 
prospective sensitivity at the end of Run II (dotted line) are also indicated (see text). 



uncertainties discussed in Sect. 01 As before, we first display our results as functions of mi/2 
moving along the WMAP strips with fixed values of Aq and tan/5. We perform a fit 
for the combined sensitivity of the observables M]y, sin^6'cfr, {g — 2)^, BR(6 — > 57), Mh 
and BR{h bb)/BR{h WW*). We do not include BR{B, ^ fi+fx-) into our fit. A 
measurement of this branching ratio at the LHC could be used in combination with the 
above measurements at the ILC. 

The results are shown in Fig. El for tan /5 = 10 and tan (3 = 50. The assumed future 
experimental central values of the observables have been chosen such that they correspond 
to the best-fit value of mi/2 in Fig. for each individual value of Aq. Thus, the minimum 
of the curve for each Aq in Fig. El occurs at = by construction. The comparison 
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Figure 17: The results of fits based on the prospective experimental accuracies for the 
precision observables Mw, sin^6'efT, {g — 2)^, BR(6 — > S7), Mh and Higgs branching ratios 
at the ILC are shown as functions of mi/2 in the CMSSM parameter space with the current 
WMAP constraints for tan/? = 10 (upper plot) and tan/5 = 50 (lower plot). For each Aq 
individually, the anticipated future experimental central values are chosen according to the 
present best-fit point. 
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of the prospective accuracies at the ILC, Fig. [TTl with the present situation, Fig. [TUl shows 
a big increase in the sensitivity to indirect effects of supersymmetric particles within the 
CMSSM obeying the current WMAP constraints. For the example shown here with best-fit 
values around mi/2 = 300 GeV (upper plot, tan/5 = 10), it is possible to constrain particle 
masses within about ±10% at the 95% CL. from the comparison of the precision data with 
the theory predictions. We find a slightly higher sensitivity for Aq < than for positive 
Aq values. For the examples with best-fit values of mi/2 in excess of 500 GeV (lower plot, 
tan/3 = 50) the constraints obtained from the fit are weaker but still very significant. 

5.2 Scan of the CMSSM parameter space 

We now investigate the combined sensitivity of the precision observables Mw, sin^ ^^efr, {g — 
2)^, BR(6 S7), Mh and BR{h hh)/BR{h WW*) in the (mi/2,Ao) plane of the 
CMSSM assuming ILC accuracies. Fig. fTHl shows the fit results for tan/3 = 10, whilst 
Fig. Uni shows the tan/3 = 50 case. 

In each figure we show two plots, where the WMAP-allowed region and the best-fit point 
according to the current situation (see Fig. [T3j) are indicated. In both plots two further 
hypothetical future 'best-fit' points have been chosen for illustration. For all the 'best-fit' 
points, the assumed central experimental values of the observables have been chosen such 
that they precisely coincide with the 'best-fit' points'^. The coloured regions correspond to 
the 68% and 90% C.L. regions around each of the 'best-fit' points according to the ILC 
accuracies. 

The comparison of Figs. ^1^1 with the result of the current fit. Fig. ^[ shows that the 
ILC experimental precision will lead to a drastic improvement in the sensitivity to mi/2 and 
Aq when comparing precision data with the CMSSM predictions. For the best-fit values of 
the current fits for tan/3 = 10 and tan/3 = 50, the ILC precision would allow one to narrow 
down the allowed CMSSM parameter space to very small regions in the (mi/2,^0) plane. 
The comparison of these indirect predictions for mi/2 and Aq with the information from the 
direct detection of supersymmetric particles would provide a stringent test of the CMSSM 
framework at the loop level. A discrepancy could indicate that supersymmetry is realised in 
a more complicated way than is assumed in the CMSSM. 

Because of the decoupling property of supersymmetric theories, the indirect constraints 

become weaker for increasing mi/2- The additional hypothetical 'best-fit' points shown in 

^We have checked exphcitly that assuming future experimental values of the observables with values 
distributed statistically around the present 'best-fit' points with the estimated future errors does not degrade 
significantly the qualities of the fits. 
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Figure 18: The results of a fit based on the prospective experimental accuracies for the 
precision observables Mw, sin^^es, {g — 2)^, BR(6 — > 37), Mh and Higgs branching ratios 
at the ILC are shown in the (mi/2, ^o) plane of the CMSSM with WMAP constraints for 
tan f5 = 10. In both plots the WMAP-allowed region and the best-fit point according to the 
current situation (see Fig. IT^ are indicated. In both plots two further hypothetical future 
'best-fit' values have been chosen for illustration. The coloured regions correspond to the 
68% and 90% CL. regions according to the ILC accuracies. 
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Figure 19: The results of a fit based on the prospective experimental accuracies for the 
precision observables Mw, sin^^es, {g — 2)^, BR(6 — > 37), Mh and Higgs branching ratios 
at the ILC are shown in the (mi/2, ^o) plane of the CMSSM with WMAP constraints for 
tan/? = 50. In both plots the WMAP-allowed region and the best-fit point for tan/3 = 50 
according to the current situation (see Fig. IT^ are indicated. In both plots two further 
hypothetical future 'best-fit' values have been chosen for illustration. The coloured regions 
correspond to the 68% and 90% CL. regions according to the ILC accuracies. 
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Figs. ^1 illustrate the indirect sensitivity to the CMSSM parameters in scenarios where 
the precision observables prefer larger values of mi/2- 

For tan/? = 10, we have investigated hypothetical 'best-fit' values for mi/2 of 500 GeV, 
700 GeV (for Aq > and Aq < 0) and 900 GeV. For mi/2 = 500 GeV, the 90% C.L. 
region in the (mi/2,^0) plane is significantly larger than for the current best-fit value of 
mi/2 ~ 300 GeV, but interesting limits can still be set on both mi/2 and Aq. For mi/2 = 
700 GeV and mi/2 = 900 GeV, the 90% C.L. region extends up to the boundary of the 
WMAP- allowed parameter space for mi/2- Even for these large values of mi/2, however, the 
precision observables (in particular the observables in the Higgs sector) still allow one to 
constrain Aq. 

For tan j3 = 50, where the WMAP-allowed region extends up to much higher values of 
^1/2 ^, "we find that for a 'best-fit' value of mi/2 as large as 1 TeV, which would lie close to 
the LHC limit and beyond the direct-detection reach of the ILC, the precision data would 
still allow one to establish an upper bound on mi/2 within the WMAP-allowed region. Thus, 
this indirect sensitivity to mi/2 could give important hints for supersymmetry searches at 
higher-energy colliders. For 'best-fit' values of mi/2 in excess of 1.5 TeV, on the other hand, 
the indirect effects of heavy sparticles become so small that they are difficult to resolve even 
with ILC accuracies. 

6 Conclusions 

We have investigated the sensitivity of precision observables, now and at the ILC, to indirect 
effects of supersymmetry within the CMSSM. We have taken into account the constraints 
from WMAP and other astrophysical and cosmological data which effectively reduces the 
dimensionality of the CMSSM parameter space. 

We have performed a analysis based on the present experimental results of the observ- 
ables Mw, sin^ ^eflf, {g — 2)^ and BR(6 — > 57) for two values of tan j3, taking into account the 
current theoretical uncertainties. For tan/5 = 10, we find that the CMSSM provides a very 
good description of the data. A clear preference can be seen for relatively small values of 
mi/2, with a best-fit value of about 300 GeV and Aq ^ —mi/2- This result can be understood 
from the separate analyses of each of the observables, each of which is well described by the 
CMSSM prediction for mi/2 ~ 300 GeV. At the 90% C.L., we find an upper bound on 
mi/2 of about 600 GeV. The supersymmetric particle spectrum corresponding to the best-fit 
region contains relatively light states. There is a possibility that some sparticles might be 
^We notice again the sparsely-populated 'voids' due to our coarse sampling procedure. 
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detectable at the Tevatron collider, and many should be detectable at the LHC [83] and the 
ILC [2] , allowing a detailed determination of their properties [93] . 

For tan /3 = 50, the quality of the fit is worse than for the case with tan P — 10. While 
Mw and sin^ ^efr prefer small values of mi/2 also for tan /3 = 50, {g — 2)^ and BR(6 — > 57) 
are better described in this case by larger mi/2 values. The indirect constraints on mi/2 are 
therefore less pronounced for tan P — 50. The best-fit value is obtained for mi/2 ~ 500 GeV 
and negative Aq. The best-fit values for the LSP mass and the lighter stau are still below 
about 250 GeV, while the preferred mass values of the heavier neutrahnos, the charginos and 
the other sleptons are in the region of 500 GeV. The 90% C.L. regions of these masses extend 
beyond 1 TeV, but would be kinematically accessible at a multi-TeV hnear coUider [94]. 
Coloured particles, such as the stops and sbottoms and the gluino are likely to have masses 
within the reach of the LHC. However, at the 90% C.L. also masses beyond ~ 3 TeV are 
possible. Heavy Higgs bosons might also be accessible at the LHC in the case of large tan /3. 

We have investigated the implications of our fit results for the prospects for detecting 
a signal for BR(S5 — > For both tan/3 = 10 and tan/3 = 50, we find that the 90% 

C.L. region for mi/2 and leads to predicted values of BR(Ss — > jjf^jJT) that are below our 
95% C.L. estimate of the Tevatron sensitivity at the end of Run II. With a more optimistic 
estimate, the Tevatron could probe a part of the parameter region for tan /3 = 50 at the 
90% C.L. It seems more hkely, however, that detection of this process would have to await 
LHC data. 

In the second part of our analysis, we have investigated the future sensitivities of the pre- 
cision observables to indirect effects of supersymmetry, assuming the experimental accuracies 
achievable at the ILC with a low-energy option running at the Z resonance and the WW 
threshold and estimating the future theoretical uncertainties. As further precision observ- 
ables besides the ones discussed for the present situation, we have included the mass of the 
hghtest CT'-even Higgs boson and the ratio of branching ratios BR(/i — > hb) /BR(/i — > 1^1^*). 
We have chosen several points in the (mi/2,^0) plane of the CMSSM with the current 
WMAP constraints as examples for 'best-fit' values, adjusting the assumed future exper- 
imental central values of the precision observables to coincide with the predictions of the 
'best-fit' values. With the prospective ILC accuracies, the sensitivity to indirect effects of 
supersymmetry improves very significantly compared to the present situation. We find that 
for assumed 'best-fit' values of mi/2 ^ 500 GeV the precision observables allow one to con- 
strain tightly mi/2 and ^40. Comparing these indirect predictions with the results from the 
direct observation of supersymmetric particles will allow a stringent consistency test of the 
model at the loop level. 
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Because of the decoupling property of supersymmetric theories, the indirect constraints 
become weaker for larger mi/2. Nevertheless, useful limits on mi/2 and Aq can be obtained for 
'best-fit' values of mi/2 as high as 1 TeV. Thus, the indirect sensitivity from the measurement 
of precision observables at the ILC may even exceed the direct search reach of the LHC and 
ILC. 

Whilst this analysis has been restricted to the CMSSM, similar conclusions are expected 
to apply if the assumption of universal soft supersymmetry-breaking scalar masses is relaxed 
for the Higgs bosons, at least for values of fi and not greatly different from those in the 
CMSSM. The impact of the dark-matter constraint may well be rather different if universal- 
ity between the soft supersymmetry-breaking squark and slepton masses is also relaxed, but 
we expect that the indication found here for relatively light sparticle masses would be main- 
tained. The investigation of these issues requires a more detailed study of models beyond 
the CMSSM, which is in preparation. 
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